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Abstract 



We discuss the dynamics of quarks within a Vlasov approach. We use an 
interquark [qq) potential consistent with the indications of Lattice QCD cal- 
culations and containing a Coulomb term, a confining part and a spin depen- 
dent term. Hadrons masses are shown to arise from the interplay of these 
three terms plus the Fermi motion and the finite masses of the quarks. The 
approach gives a lower and an upper bound for hadrons. The theoretical 
predictions are shown to be in fairly good agreement with the experimental 
data. 
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Very recently the approval and construction of new machines capable of accelerating 
light and heavy ions at 100 GeV/nucleon and more has stimulated the interest for the 
search of a quark gluon plasma [|I|. Methods developed to describe the dynamics of two 
colliding heavy ions at low energies can be very useful in this new area when modified to 
take into account the intrinsic structure of the hadrons. As an example we discuss in this 
paper an approach based on the Vlasov equation (VE). This equation can be obtained from 
the quantum Bogoliubov-Born-Green-Kirkwood-Yvon(BBGKY) hierarchy by means of the 
Wigner transform in the limit h ^ It can include quantum effects such as the Pauli 

principle and can be easily extended to relativistic dynamics [0,0. In this paper we will 
apply the VE to the dynamics of the constituents quarks in hadrons. The quark masses 
can be quite small as in the case of (u,d,s) quarks, (see table I for the values used in this 
paper), so that a relativistic treatment is necessary. As far as the interaction is concerned 
we will start from a two body potential consistent with Lattice Quantum Chromo Dynamics 
(LQCD) calculations [|T],0, which displays both a Coulomb U{r — >■ 0) oc 1/r behaviour and a 
confining part U{r ^ oo) oo. In particular, we will use the Richardson's potential which 
depends on the scale parameter A From such a starting point we calculate the time 
evolution of the quarks and the masses of the hadrons. The general experimental features 
are quite reasonably reproduced, even though the radii of light hadrons are underestimated. 
We notice that a similar approach has been proposed in 0, but a different interaction was 
used. 

We briefly recall some general features of the VE, interested people should look the (not 
complete) list of references for more details on this equation The VE gives the time 

evolution of the one body distribution function fq^ (r, p, t) in phase-space: 

d,u^ + ^ ■ V./,, - v,u ■ Vp/,, = (1) 

where U(r) is the potential (discussed below) which governs the quarks dynamics, 
E = \/'p^ + VP? is the energy and m is the quark mass. The underscripts indicate that 
the distribution function depends on the flavor (q) and color (c) of the quarks. Notice that 
the kinetic part is properly relativistically treated, while the potential term is non relativis- 
tic. A relativistic extension of the Richardson's potential is discussed in but will not be 
included here. 

Numerically the VE equation is solved by writing the one body distribution function as: 

1 ^ 

fqcir,P,t) = — Y.9r{r - ri{t))gp{p - p,{t)) (2) 

where the gr and gp are sharply peaked distributions (such as delta functions, gaussian 
or other simple functions), that we shall treat as delta functions. = Qutp is the number 
of such terms, Q = q + q is the total number of quarks and antiquarks (for a meson Q=2). 
Actually, N is much larger than the total quark number Q, so that we can say that each 
quark is represented by nxp terms called test particles(tp). The rigorous mean field limit 
can be obtained for utp ^ oo where the calculations are of course numerically impossible, 
even though the numerical results converge rather quickly. Inserting eq.( |^) in the Vlasov 
equation gives the Hamilton equations of motions for the t.p. 0]: 
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Pi 

Pi = -Vr,U, (3) 

for i = 1....N. The total number of tp used in this work ranges from 5000 to 50000 with 
no appreciable change in the results. The equations of motion (3) are solved by using a 
0((5t^) Adams-Bashfort method 

Let us now specify the gg-potential. In agreement to LQCD calculations we have 
for mesons(^ = 1): 

U(r) = A(Ar - ^^-V^) + — — — 5(r) (4) 



where 



n/ is the number of quark flavors involved and the parameter A has been fixed to re- 
produce the masses of heavy cc and bb systems in 0. In eq.(4) we have added to the 
Richardson's potential the chromomagnetic term, very important to explain the masses of 
different resonances for light quarks. In this work the expectation value of < aga^ > is 
used depending on the relative spin orientations of the constituent quarks . For instance for 
the pion this term is equal to -3 while for a p meson it is equal to +1 [Q. The 6 function 

is approximated to a gaussian i.e. we make the replacement 6{r) — > ^2-Ka^)3/^ ^~^ ^ ^^'^ 
we fixed a = 0.5/m. The results depend on the ratio of the average value of the strong 
coupling constant ag to the variance of the gaussian. Such ratio was fixed to fit the mass 
difference between n and p mesons (see below). Clearly the chromomagnetic interaction 
becomes unimportant for very heavy quarks. The potential term eq.(4) acts between two 
quarks, since we are describing quarks as a swarm of tp, we normalize the potential by a 
factor l/ritp 0. 

The initial conditions are given by randomly distributing the tp in a sphere of radius 
r in coordinate space and in momentum space. P/ is the Fermi momentum estimated 
in a simple Fermi gas model by imposing that a cell in phase space of size h = 2tt can 
accommodate at most two identical quarks of different spins. A simple estimate gives the 
following relation between the quarks density Ug and the Fermi momentum: 



n. 
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V) (6) 



an analogous formula can be derived for q 0. The degeneracy number = -n-c x x n^, 
where Uc is the number of colors and Ug is the number of spins . For quarks and antiquarks 
3 different colors are used red, green and blue (r,g,b) From the above equation we see 
that the Fermi momentum for quarks distributed in a sphere of radius 0.5 fm is of the order 
of 0.5 GeV/c. Thus relativistic effects become important for quark masses less than 1 GeV. 
For instance for the vr case discussed below, if we calculate the total energy of the system 
relativistically we obtain 0.14 GeV, while using the nonrelativistic limit we obtain about 1 
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GeV! We stress that since the system is properly antisymmetrized at time t=Ofm/c, it will 
remain so at all times since the VE conserves the volume in phase-space [§]. 

The masses of the hadrons are determined by finding the minumim total energy of the 
system as a function of the initial radius r. For each initial radius r a Fermi momentum is 
deduced from eq.(6) and this gives (after adding the potential term) a total energy of the 
system which is of course constant in time. Changing the initial radius changes the total 
energy, and it has a minimum at r = f. Thus it is the interplay among the Fermi motion, 
the potential term and the quark masses which determines the masses of the hadrons. 

In order to test our approach, we have first studied heavy quark systems using the same 
values for the c and b masses, nf = 3 and A = 0.398Gel^ as in ref. and neglecting the 
chromomagnetic term. We obtained the minimum mass values of 2.9 and 9.5 GeV for the 
cc and bb systems, to be contrasted with 3.1 GeV and 9.4 GeV obtained in 0. The good 
agreement to the quantum calculation of ^ suggests that the VE is quite well justified and 
in particular the mean field approximation is good. This is probably due to the fact that 
the interparticle potential is the result of many gluons exchanges and that the hadrons are 
made of a combination of colored quarks and in this sense the mean field approximation is 
reasonable. Also, it is important to stress that in a quantum calculation the bare two body 
potential is folded with smooth qq wave-functions. The smoothness of the wave functions is 
simulated in the VE with the use of (a large number of) tp. 

In order to study the limits and merits of our approach we have extended the calculations 
to lighter quarks. For such systems the chromomagnetic term is important to reproduce 
the experimental values of the masses. Using the same value of the scale parameter A = 
0.398(761^ as in P and fixing the (u,d) quark masses and the strong coupling constant as, 
we can easily reproduce the pion mass, but not the p mass at the same time. Thus we have 
readjusted the values of the scale constant and the quark masses to reproduce the data. We 
found a good fit for mesons by using A = 0.250GeV, as = 0.225 |Q and the quark mass 
values given in Table I. We notice that the parameters and the heavy quark mass values 
are in good agreement with currently accepted ones The masses of (u,d) quarks in 

table I is somewhat smaller than the 300 MeV used in many potential models 0]. This is 
due to the fact that these models are non relativistic but, as we stressed above, because 
of the Fermi motion, relativistic effects are quite important. In the relativistic approach of 
0, where the Richardson's potential was used as well, the (u,d,s) quark masses have values 
comparable to ours. The small differences between our results and are most probably 
due to their relativistic generalization of the Richardson's potential. 

The parameters entering our model are essentially fixed on some meson masses. As a con- 
sistency check we extended the calculations to the baryons case with the usual modifications 
of a factor 1/2 to the potential eq.(4), as suggested by LQCD considerations Q. 

From the knowledge of the distribution function at each time step, we can easily calculate 
the density and the potential of the system at one time step. An average over time of these 
quantities can be performed as well. In the calculations the time averages were performed 
over a time interval up to 100 fm/c. In figure (1), we plot the density (left column) as a 
function of the distance r from the center of the system and obtained as average over the 
one body distribution function at one time step(square symbols), and over time as well (full 
line). The top panel correspond to a total energy of a (u,d) system of 0.140 GeV and a root 
mean square (RMS) radius (averaged over time) of 0.33 fm, i.e. the pion. We stress that 
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this is the minimum for the total energy obtained when the total quarks spin add to S=0, 
i.e. a pseudoscalar mesons [^. The middle panel corresponds to a total energy of 0.85 GeV 
and a RMS of .53 fm for a (uud) system, i.e. a nucleon. A similar calculation but for S=3/2 
(A) gives a mass of 1.3 GeV and a RMS of 0.75 fm (bottom panel in figure 1). In all cases, 
time and tp averages are identical which implies that the systems are in equilibrium. The 
densities obtained averaging over the distribution function at one time step, do not extend 
to very low values because of the finite number of tp. The rms radius of the vr and n are 
smaller than data 0], however the densities extend much further than the rms radius, and 
they fall off exponentially similarly to experiments. For increasing quark masses the rms 
becomes more reasonable as compared to data or other calculations. 
The average potential can be easily estimated as 

U{r) = —Y^U{r,r,) (7) 

rifp j 

The average potential displays some interesting features, see fig. 1 (right column). First, 
all the potentials go to zero for all systems (even though it is explicitly shown in the figure 
for the A-case only), and for large distances in contrast to the bare potential eq.(4) that 
diverges linearly, just because the densities go to zero. The average potential for vr and n is 
an attractive pocket, and the confining term gives some contribution at large distances. For 
the A case, the chromomagnetic term is repulsive thus the RMS of the system is larger than 
n and the confining term only is responsible for keeping the system bound. Looking at the 
A-density we see that the repulsion at small distances gives a smaller and flatter density as 
compared to the n case. 

In figure (2) we display the calculated (open symbols) mass of the resonances vs. the 
sum of the quark masses (cfr. table I), for mesons (top) and baryons (bottom). The circle 
symbols refer to attractive, while the squares to repulsive chromomagnetic term in eq.(4) 
The corresponding experimental data [1^ are given by the full symbols. The overall 



agreement is quite good in all cases and some predictions for resonances not yet observed 
are also given. 

In order to understand how other higher resonances appear for fixed quark types, we 
define a dipole operator for mesons, analogous to the nuclear case PJTT[]: 

D{t) = Y.vAt) -Y.vAt) (8) 

where the sum is extended over all the tp for q and g, the choice of the z-axis is of course 
arbitrary. Resonances are better seen by defining the dipole strength function S(E), where 
E is the energy: 

S{E) = \Fm^ (9) 

with the Fourier transform 

F{E) = J dt[e'^^'^^^D{t)] (10) 

In figure 3 we plot the function S(E) vs E for the vr (ud-quarks) case. Strong resonances 
extend to about 1. GeV. This implies that we can have mesons up to 1. GeV above the 
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ground state mass. In the case of the heavier bb quarks, resonances extend to about 0.5 
GeV above the ground state. Experimentally resonances are seen for instance for the bb 
case to about 2 GeV above the smallest resonance [|10|. In our model we find that to obtain 
resonances at higher energies we need to increase the strength of the confining term. Thus 
the measured upper values of the hadron masses can give a constraint on the value of the 
string tension. 

The results discussed above prove that the VE is suitable to describe the quarks dynamics 
in hadrons despite the simple form for the two body potential. The approach works rather 
well for heavy quarks, in that it is able to reproduce the masses and radii of the heavy 
hadrons (as compared to data or other calculations 0,^) and also the masses of light hadrons. 
The radii of light quarks systems are underestimated which could be a hint for relativistic 
corrections to the potential term as discussed for instance in [|^. 

The method is quite easy to implement and the numeric is rather well under control. The 
equation of state of quark matter can be calculated within the same formalism. Future work 
will be also to implement a collision term which should help to understand the dynamics of 
colliding hadrons. 
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Table I 



Quark masses used in this work. 



QuarkMass 


GeV 


u 


0.13 


d 


0.13 


s 


0.35 


c 


1.45 


b 


4.8 


t 


180. 



8 



FIGURES 




FIG. 1. Average density (left column) for a -7r(top), ?7-(middle) and a A (bottom) hadrons. The 
average is over tp (squares) and over tp and time (full line). Similarly for the average potential 
(right column). 
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FIG. 2. Meson masses (top) vs. masses of the qq pair. The symbols refer to the observed 
(full symbols) and calculated (open symbols) masses. The chxles refer to the pseudoscalar and the 
squares to vector mesons masses. Similarly for baryons (bottom). Data are taken from [9]. 
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FIG. 3. Strength function (in arbitrary units) versus energy for the tt case. 
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